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Abstract 

We find the equation of state p, p oc T 6 , which gives the value of the sound 
velocity c 2 = 0.20, in agreement with the "realistic" equation of state for hot 
hadronic matter suggested by Shuryak, in the framework of a covariant rela- 
tivistic statistical mechanics of an event-anti-event system with small chemical 
and mass potentials. The relativistic mass distribution for such a system is 
obtained and shown to be a good candidate for fitting hadronic resonances, 
in agreement with the phenomenological models of Hagedorn, Shuryak, et al. 
This distribution provides a correction to the value of specific heat 3/2, of the 
order of 5.5%, at low temperatures. 
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1 Introduction 

One of the main goals of experiments with high-energy nuclear collisions is to produce 
and to study hadronic matter, in particular, trying to reach conditions at which the 
phase transitions into the quark-gluon plasma phase can take place [l]-[5]. The physics 
of hot hadronic matter has not been studied much, although such matter is already 
produced in current experiments in high-energy physics. 

The experimental data on multiple hadron production obtained in recent years are 
in agreement with the main consequences of the theory formulated by Landau || over 
40 years ago. However, with a sufficiently quantitative approach, it becomes necessary 
to consider a number of physical effects which bring about certain modifications 
of the results obtained in the fundamental work [6]. For instance, the solution of 
the equations of motion obtained by Landau differs from a numerical calculation of 
Melekhin 0, as a result of an inaccurate estimate. A more accurate analytic solution 
was given in [|J. 

The equation of state in Landau's work is taken to be p = p/3 (where p is the 
pressure and p is the energy density), corresponding to an ultra-relativistic gas. How- 
ever, the lab energies s = {P1+P2) 2 ~ 3x 10 3 GeV 2 correspond to initial temperatures 
T ~ 1 GeV. In this temperature range, the interaction of the hadrons is strong and 
has mainly a resonant character, the masses of the resonances being comparable with 
the temperature. Thus, hadronic matter under these conditions is neither an ideal 
nor an ultra-relativistic gas. 

Corrections in the equation of state due to the interaction of the hadrons have been 
discussed in the literature for the last three decades [1],[4],[8]-[14]. These considera- 
tions were based mostly on a phenomenological model in which the Landau theory is 
applied to all particles except the leading ones, i.e., the fragments of the initial parti- 
cles, whose characteristics are taken directly from experimental data. The framework 
for the latter considerations is the QCD phase transition of hadronic matter into the 
quark-gluon plasma. 

Hot hadronic matter, in which volumes per particle are a few cubic fermis, is 
certainly made out of individual hadrons. It is clear that, at low temperatures T << 
m n , one has a very rare (and therefore ideal) gas of the lightest hadrons, the pions. 
As the temperature is raised and the gas becomes more dense, one should take into 
account interactions among the particles. This was done using the following three 
approaches: (i) the low-temperature expansion, (ii) the resonance gas, (iii) the quasi- 
particle gas. 



2 



The first approach for the pion gas is based on the Weinberg theory of pion inter- 
actions [1511) which uses the non-linear Lagrangian containing all processes quadratic 
in pion momentum. Its first application to calculation of the thermodynamic param- 
eters of the pion gas was done in ref. [1]. One of the important consequences of this 
work was that, after isospin averaging, all corrections quadratic in momenta cancel 
each other, and corrections proportional to the pion mass (Weinberg tttt scattering 
lengths) are nearly compensated, producing negligible corrections at the 1% level. 
Corrections of the second order in the Weinberg Lagrangian were also estimated in 
[1], and a much more systematic study of the problem including quartic terms in the 
mesonic Lagrangian were made in ref. |I6 . Without going into discussion of these 



works we remark that the applicability of this approach is limited by temperatures 
T < 100 MeV, for which typical collision energies are significantly below resonance 
.However, such T are lower than even the lowest temperature available in experiments, 
because the so-called break-up temperatures are typically T ~ 120 — 150 MeV. 

The idea of resonance gas was first suggested by Belenky and Landau [|I7[], as early 
as 1956. They used the Beth-Uhlenbeck formula [|Tj|, well known in the theory of the 
non-ideal gases ||19|| , which relates the second virial coefficient with scattering phase 
shifts Si(p). In particular, in this approximation the "internal" part of the statistical 
sum can be written as a sum over bound states plus the scattering part[|: 

Z** = E z- EJT + - E / dp^e~^ T . (1.1) 

Using Eq. (1.1) one can prove that a narrow and elastic resonance is analogous to 
adding one extra physical state to the system. 

This observation was later used by many authors, in particular, by Hagedorn 
0, who has generalized this statement to the idea that all hadrons (both stable 
and resonances) should be treated as real degrees of freedom of hadronic matter. The 
statistical bootstrap model Pl,fl0|l initiated by Hagedorn, which requires that the 



spectra of the resonances and of the system as a whole coincide, provides a resonance 
spectrum 

p{m) ~ m a exp(m/T ), (1.2) 

where a and To are some parameters. The statistical sum diverges when T > To, which 
indicates that the theory involves a limiting temperature (the so-called "hadronic 
boiling point" |J) to which the system can be heated. Such behavior indicates that 
there is some phase transition, and that the language used is meaningless above it. 
This phase transition, as well as thermodynamic consequences of the mass spectrum 
(1.2), were discussed in ref. [Q. The expression (1.2) with a = —5/2 is in good 
agreement with experiment in the resonance mass region m ~ 1.2 GeV M. 

lr This formula is written down for the nonrelativistic case, as it appears in textbooks on statistical 
mechanics, e.g., in u9|. 
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In 1972, Shuryak ||, in place of (1.2), used the simple power parametrization 



p(m)~m fc , (1.3) 

which also describes experiment in this mass region for k ~ 3. The energy density is 
expressed by [|J [E 2 = p 2 + m 2 ) 

P = E/ 'dmp(m) f a = ±l. (1.4) 

<j J (2vr) d exp(y) + a 

Since p(m) is a rapidly growing function, the main contribution to integrals of type 
(1.4) is given by the mass region in which exp(m/T) >> 1, and the difference in 
properties of bosons and fermions is irrelevant. Upon substitution of (1.3) into (1.4) 
one finds that the energy density and, analogously, the pressure, defined by 

P= lj2 fdmp(m) ,L ~ > ( L5 ) 

o„J J {Aix) A h, exp(^) + a 

is proportional to T k+5 . Therefore, the velocity of sound c, defined by 



c = T P > (L6) 

turns out to be a temperature independent constant, and the equation of state belongs 
to the class considered by Melekhin For this case, the expressions for energy and 



pressure are as follows, 

p = XT k+5 , p =-^T k+5 , (1.7) 
r k+4 v ' 

where A is a certain constant. Consequently, for k ~ 3 [8] 

c 2 = — — w 0.14. (1.8) 
k + 4 v ' 

This numerical value of c 2 , although less than the ideal gas value c 2 = 1/3, is obtained 
very approximately due to replacement of the sum over the resonances by the integral 
(1.4). Strange and baryonic resonances may only occur in pairs, which is not taken 
into account in derivation of (1.8). Moreover, at temperatures T ~ m, which is the 
case we are interested in, quantum field theory requires the inclusion of antiparticles 
which becomes important and therefore should be taken into account. 

In 1975, Zhirov and Shuryak |11[] used the Beth-Uhlenbeck method |18|| , which 
reduces the problem in the case of a purely resonant interaction to a simple case 
of a gas consisting of a mixture of stable and unstable particle (resonances) treated 
on an equal footing. Their calculations give the value of the velocity of sound c 2 = 
0.18 — 0.21, for the range of temperatures 0.2 — 1.0 GeV, which is appreciably smaller 
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than than the asymptotic value c 2 = 1/3 and larger than the estimate c 2 ~ 0.14 of 
ref. [1. 

In his book ||] of 1988, Shuryak proposed the "realistic" equation of state (referred 
to also in [1]) 

p ~ (20GeV^ 2 ) T 6 , (1.9) 

describing the behavior of hadronic matter in the temperature range 0.2GeV ~ T ~ 
l.OGeV (which gives c 2 = 0.20), in good agreement with the relevant experimental 
data. 

A similar equation was obtained in refs. ||22||,|[23|| within the concept of hadronic- 
plasma phase transitions, for a kind of chiral bag model for the quark-gluon plasma 
cluster. 

In ref. Shuryak, in order to describe the behavior of the energy density with 
temperature p oc T 6 , proposed to consider hot hadronic matter as a gas of quasi- 
particles, which have quantum numbers of the original mesons, but with dispersion 
relations modified by the interaction with matter (similarly to Landau's idea of "ro- 
tons" , to explain the growth of the energy density of liquid 4 He with temperature 
more rapidly than T 4 ). Since in such an approach hadrons are included as physical 
degrees of freedom of hot hadronic matter, they are assumed not to be absorbed too 



strongly, i.e., to be good quasiparticles. In ref. [12] an attempt was made to guess 



what these dispersion relations should be, in order to obtain the expected behav- 



ior of the energy density with temperature. In |14| these dispersion relations were 
calculated in explicit form. 

In this paper we show that the same results can be obtained within the framework 
of a manifestly covariant relativistic statistical mechanics [24]- [27]. We shall review 
this framework briefly in the next section. First a possible correlation of the manifestly 
covariant theory with the hadronic spectrum of the Hagedorn-Frautschi form was 
discussed by Miller and Suhonen [28] . They remarked that an approach based on the 
grand canonical distribution function of ref. [24] gives results similar to the ones for 
the resonance gas in the case ~ 0, where fix is the additional mass potential of 
the ensemble [24] . In the present paper we study a free relativistic gas in the limiting 
casef] — > oo (where M is an intrinsic scale parameter for the motion in space- 
time [29]), taking account of antiparticles. For such a particle-antiparticle system 
with n ~ (i.e., zero net particle charge) we find the equation of state p, p cx T 6 , 
which gives the value of the velocity of sound c 2 = 0.20. We obtain the relativistic 
mass distribution for this system and show that this distribution is a good candidate 
for fitting the hadronic resonances, in agreement with the theory of Hagedorn [9], [10]. 



2 For finite temperature this case is realized by \ik —> 0. 
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2 Relativistic TV-body system 



In the framework of a manifestly covariant relativistic statistical mechanics, the dy- 
namical evolution of a system of N particles, for the classical case, is governed by 
equations of motion that are of the form of Hamilton equations for the motion of N 
events which generate the space-time trajectories (particle world lines) as functions 
of a continuous Poincare- invariant parameter r, called the historical time [29], [30]. 
These events are characterized by their positions = (t, q) and energy-momenta 

= (E,p) in an 8iV-dimensional phase-space. For the quantum case, the system 
is characterized by the wave function ip T (qi, q2, ■ ■ ■ , <Zjv) £ L 2 (R 4N ), with the measure 
d A qid A q2 ■ ■ ■ d A q^ = d m q, (<& = gf ; /i — 0, 1, 2, 3; i — 1, 2, . . . , N), describing the dis- 
tribution of events, which evolves with a generalized Schrodinger equation [29]. The 
collection of events (called "concatenation" [31]) along each world line corresponds 
to a particle, and hence, the evolution of the state of the iV-event system describes, 
a posteriori, the history in space and time of an iV-particle system. 

For a system of N interacting events (and hence, particles) one takes [31] 

K = Y<YM+n<luq2,...,qN), (2.1) 

i 

where M is a given fixed parameter (an intrinsic property of the particles), with the 
dimension of mass, taken to be the same for all the particles of the system. The 
Hamilton equations are 

dq? _dK _pt 
dr dp ifl M' 

d$ OK dV 
dr dq ifX dq ifX ' 
In the quantum theory, the generalized Schrodinger equation 



(2.2) 



d 

i-jf-^AQi, 52, • • • , Qn) = Kip T (q u q 2 , . . . , q N ) (2.3) 
describes the evolution of the iV-body wave function ip T (qi, Q2, ■ ■ ■ , Qn)- 

2.1 Ideal relativistic identical system 

To describe an ideal gas of events obeying Bose-Einstein/Fermi-Dirac statistics in the 
grand canonical ensemble, we use the expression for the number of events found in 
[24], 

N = V n ^ = Ve( E -^«^y T Tl ^ 
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where px is the additional mass potential [24], and m 2 = —k 2 = —k^k^. Replacing the 
sum over k^ by an integral, one obtains for the density of events per unit space-time 
volume n = N/V^ [27], 

1 r°° m 3 dm sinh 2 3 dB , 
n=— 2 M ■ 2.5 

All 3 Jo g(m cosh P-pl+^k )/T -p i 

The integration results in the following expression for n [27]: 

1 M 2 00 pt sM 

■ T 2 E(±l) s+1 ^r^(2,2;-^-), (2.6) 



(2tt) 3 /4 ^ v 7 s 2 v ' '2^T" 



where ^/(a,6;z) is the confluent hypergeometric function 

The formula (2.4) should be considered as the analytic continuation of (2.6). As 
a power series, (2.6) diverges for |e M//T | > 1, i.e., for p > 0, while (2.4) is always valid. 

In the limit » 1 (or T « -^), Eq. (2.6) reduces to [27] 

n = ±^T 4 L* 4 (e^), (2.7) 

fc 

where Li n (z) = J2kLi frr is the polylogarithm [33]. One can also obtain in this limit, 

P = ±^7*Li 6 (± e *) = 5p, 
iVo = ±J^T 5 L, 5 (±e^), (2.8) 



7T 



where p is the pressure and p is the density of energy of the particle gas, and Nq is the 
density of particles per unit space volume; Tav is the average passage interval in r for 
the events which pass through a small (typical) four-volume AV in the neighborhood 
of the i? 4 -point [25]. 



3 Introduction of antiparticles 

In establishing the theory presented in this paper we assumed that the total number 
of events and, therefore, particles, is a conserved quantity, so it makes sense to talk of 
a box of iV particles. This can no longer be true at high temperatures [34]; it is well 
known that at temperatures T ~ m quantum field theory requires the inclusion of 
particle-antiparticle pair production, which becomes important and therefore should 
be taken into account. If N is the number of antiparticles, then N and N by them- 
selves are not conserved but N — N is. Therefore, the high-temperature limit of (2.4) 
is not relevant in realistic physical systems. 
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The introduction of antiparticles into the theory in a systematic way was made by 
Haber and Weldon [34] within the framework of the usual on-shell theory. They con- 
sidered an ideal Bose gas with a conserved quantum number (referred to as "charge" ) 
Q, which corresponds to a quantum mechanical particle number operator commuting 
with the Hamiltonian H. All thermodynamic quantities may be then obtained from 
the grand partition function Tr {exp[— (H — Q)/ksT]} considered as a function of 
T, V, and fi []35f| . The formula for the conserved net charge readsf] [34] 



k 



e {E k -^)/k B T _ I e (E k +n)/k B T _ I 



(3.1) 



In such a formulation a boson-antiboson system is described by only one chemical 
potential /i; the sign of /i indicates whether particles outnumber antiparticles or vice 
versa. For the case of fermions one has simply to change the sign in the denominators 
ofEq. (3.1). 

The introduction of antiparticles into the theory we are discussing here, as the 
events in the CPT-conjugate state, having nonnegative energy, leads, by application 
of the arguments of Haber and Weldon, to a change in sign of fi in the distribution 
function for antiparticles. A full consideration of particle-antiparticle pair production 
within the off-shell framework should include anti-events as well, i.e., events having 
the opposite sign of \ik- 



The full theory of anti-events is constructed in [Q. Upon r reversal combined 
by charge conjugation, they are included in statistical mechanics as the usual events 
(with positive M), having the opposite sign of fix in the distribution function p6|,[37]. 
The following relation which generalizes (2.4) is found as the analogue of the formula 
(3.1): 

N = N E + N E -N E -N E 

= N(fi, /i K ) + N(-/jl, -/i K ) - N(-/i, fi K ) - N(n, -ji K ) 

= E[n fc M(/i,/ix) + n kf >(-n, -Hk) - "^(-ftM - nkn(ji, -^k)} 



E 

kv 



1 1 

+ 



e (E-^ K ^)/k B T -p x e (E+v-n K %M)/k B T -p l 
1 1 



(3.2) 



where N is the conserved net event "charge" , Ne and N E are the numbers of events 
with E > 0, E < 0, respectively, and Ne, N e are the corresponding numbers of 
anti-events. 



3 One uses the standard recipe according to which all additive thermodynamic quantities are 
reversed for antiparticles. 
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Replacing summation by integration, we now obtain a similar formula for the 
event number densities: 



n 



47T 3 JO 



m 3 dm sinh 2 (3d(3 



1 



+ 



1 



e (mcosh/3-/i+^ K 2L_)/T -p ^ g (m cosh/3+^-^t /f ||j )/T -j- j 
1 1 



e (mcosh/3+ju+Aiisr^f)/T -p ^ e (mcosh / 3-/x-/i K f| i -)/T -j- j 

(3.3) 

The event-anti-event system possessing Bose-Einstein statistics was treated in ref. 
[37]. The requirement that the four n^'s in Eq. (3.2) be positive definite leads to 
the bounds on the mass spectrum [37] 



M 



1 



and to the relation 



M 



< m < 



\Wk\ < 



M 
\Hk\ 

M 



M 



The mass region (3.4) for small /i is well approximated by 

2M 

\fjL\ < m < 



(3.4) 



(3.5) 



(3.6) 



Therefore, for \fi\ ~ and \hk\ —>■ we have almost the whole range of m, (0, oo). 
Since there is no requirement on the n^'s in the fermionic case, we take for this case 
< m < oo. Thus, both cases can be treated simultaneously, with the whole range 
of m. In this way we obtain, upon expansion of the denominators in Eq. (3.3) into 
power series and integration on /3, 



T r°° x % (±l) s+1 2 . sfi . sjiKmr (sm 

— I dm > m sinh — smn Ki — 

7T 3 Jo ^ s T 2MT 1 V T 



n 



s T 2MT 

from which we identify the mass distribution f(m), normalized a; 

dmf(m) = \n\. 



In (3.7), K\ is a Bessel function of the third kind (imaginary argument). 
Integration in (3.7) gives the following expression for n : 



(3.7) 



(3.8) 



1 M 2 ~ (iir+1 . an 



n 



s=l 



*(2,2; 



*(2,2; 



2fx K T y 



(3.9) 



We remark that only |n| has physical meaning as a conserved net event charge. 
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The formulas for p, p and Nq become []36 



P 



P 



o 



p(p, Pk) + p(~P, ~Pk) + p(~P, Pk) + p(p, ~Pk), 
p{p, Pk) + p(~P, ~Pk) + p{~P, Pk) + p{p, ~Pk), 
N (p, Pk) + N (fi, -jj, K ) ~ N (-[M, p, K ) - N (-fj,, -p K ) 



(3.10) 
(3.11) 
(3.12) 



where p(/z, ^k), p(p,Pk) and N (fi,fi K ) are the subdistributions defined by the cor- 
responding formulas of the event statistical mechanics [27]. 

Note that the expressions (3.12) for N (actually corresponding to the Q in Eq. 
(3.1)), which refers to particles and antiparticles, and (3.3) for n, which refers to 
events, differ in sign of the corresponding terms containing —fix- While iVo com- 
mutes with the effective Hamiltonian of a relativistic many body system in the usual 
framework, it is (3.2) which commutes with the evolution operator of the covariant 
theory [0ql . 

We also remark that, as in the usual theory, the system we are studying is de- 
scribed by only one chemical potential fi and only one mass potential the sign of 
\i indicates whether particles outnumber antiparticles or vice versa, while the sign of 
\x\xk has a similar relation to the relative number of events and anti-events, as seen 
in Eqs. (4.8). 

4 Realistic equation of state for event— anti-event 
system 

In the limit >> T, we use the asymptotic formula for z — » oo (ref. [32], p. 278, 
subsection 6.13.1) 



\ z z° 
and obtain from (3.9), retaining the terms up to the order of ( i ^-) 2 , 




) 



(4.1) 



n 




(4.2) 



(4.3) 



The corresponding formulas for p, p and N read [27] 



P 
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M 2 
4 T AV 



+48^ [Li 8 (±e$) + Li 8 (±e-$) 
T 2 {5 [Li 6 (±er) +Li 6 (±e 



7T 3 M 



+432^ [Li 8 (±e$) + Li 8 (±e-*) 



± 



M 2 
4 T Ay 



7T 3 M 

+48 



/4g 

M 2 



T 5 {[Li 5 (±e^)-Li 5 (±e-^) 
Li 7 (±eT) — Li 7 (±e~T) 



Note the remarkable relation 



AY 



v 



\Pk\ 



\n\. 



(4-4) 



(4.5) 



(4.6) 



(4.7) 



First consider the bosonic case. Using the corresponding relations for the polyloga- 
rithms (see Appendix A), and retaining the terms up to the order of (^fp) 2 , (^) 2 , 
we obtain 



n = 



V = 



47T lifiKm* 



45 M 



rpq 



5 /i 2 96w 2 p 2 K T 2 



2tt 2 T 2 



21 M 2 



1tt 3 T, 



945 M 



/ 1 + ^l/^ + 24vr 2 /i 2 f T 2 



P 
iV 



407r 3 T Aym6 



945 M 
47rT Ay 



4tt 2 T 2 
' 21 /i 2 



4tt 2 T 2 



5 M 2 J 
216n 2 p 2 K T 2 ' 
25 M 2 



45 M 



/iTM 1 + 



5 /i 2 96tt 2 /4T 2 



2tt 2 T 2 



+ 



21 M 2 



Therefore, 



c 2 = ^^l(l 
dp p 5 V 



96 /4g 
25 M 2 



(4.8) 



(4.9) 



i.e., one obtains the value of the velocity of sound c 2 = 0.20, as well as the correction 
to this value, of the order of { li ff) 2 - 

In the fermionic case, one uses the corresponding relations for the polylogarithms 
with the argument z < (Appendix A), and obtains similar formulas for n, p, p and 
No, which give 



c 2 - 1 - 



127 96 fi? K (k B T)' 



124 25 



M 2 



We see that in the case p ~ 0, px — > the equation of state for event-anti-event 
system reads p, p oc T 6 , implying the velocity of sound c 2 = 0.20, in agreement 
with the Shuryak's "realistic" equation of state. We now turn to the corresponding 
distribution of mass. 
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5 Relativistic mass distribution in event— anti-event 
system 

The distribution of mass, defined by Eq. (3.7), can be simplified in the case /i ~ 0, 
rr—jy >> 1- Since for the whole range of m and small /i the average (m) depends 

only on temperature [26], [27], ^*mt * s °^ ^ ne OT( ^ eT °f << 1- Therefore, one can 
replace the hyperbolic sine functions by their arguments and obtain in this way 

1 \nn K \™ +1 4 fsm\ 
t(m) ~ — 5-——=- > ±1 sm Ai — — . 5.1 
JK J 2vr 3 MT ^ V T J 17 

For a given T, in the region m << T we use the asymptotic formula 
(ref. f38|, p.375) 

and obtain 

/(m)~m 3 , (5.2) 

for both bosonic and fermionic cases. If we assume that the distribution (5.2) de- 
scribes the genuine spectrum of compound hadron-resonance systems, and restrict 
ourselves only to the region m « T, we shall obtain, with Shuryak's Eqs. (1.3), (1.8) 
for k = 3, 

c 2 = 0.14. 

In the region m » T, we use another asymptotic formula (ref. P8| , p. 378) 




Since in this region particles become distinguishable, we have 

f(m) ~ m 7/2 e"^. (5.3) 

In this case the equation of state tends asymptotically to the equation of state of an 
ideal Stefan-Boltzmann gas, p = p/3, providing for the velocity of sound the value 

c 2 = 0.33. 

For the transitional region m ~ T, the formula for K v (z) can be found by the method 
of steepest descents (ref. [32], Vol.2, p. 28). The corresponding expression for the 
mass distribution can be also obtained directly from Eq. (5.3), since in this case the 
power function dominates over the exponent. In either case, one has 

f(m) ~ m 7/2 , (5.4) 
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so that, if we restrict ourselves only to this region, we shall obtain, through Eqs. 
(1.3), (1.8) with k = 7/2, 

9 2 

c 2 = — sa 0.13. 
15 

If the three cases are treated on an equal footing, one has 

c 2 = -(0.14 + 0.13 + 0.33) = 0.20, (5.5) 
3 

in agreement with the equation of state (4.9) for the whole mass region. We see that 
the mass distribution (5.1) is a good candidate for fitting the hadronic resonances, in 
good agreement with the phenomenological models of Hagedorn, Shuryak, et al. 



5.1 Specific heat of event— ant i- event gas 



In conclusion we wish to calculate the specific heat of the event-anti-event gas we are 
studying. The distribution (5.1) gives the following value of the average m : 



(m) 



Jdmmf(m) _ J^dmm^^f) 



457T 



fdmf(m) f£°dmm*Ki(f) 32 
where we used the formula (ref. [39], p. 684, formula 16) 



T. 



(5.6) 



dx x^KJax) = 2^- 1 a->*- 1 r 



1 + fi + u 



r 



l + fi—U 



(5.7) 



The rest frame average of the energy (mcosh/5) can be obtained from (3. 3), (5.1) 
with the help of the relation 

sinh 2 (3 cosh (3 = - (cosh 3/5 — cosh (3) 



and the formula (ref. [39], p. 358, formula 4) 



dx cosh vx e 



-a cosh x 



KJa) 



(5i 



as follows, 

(E) 



m 3 dm sinh 2 (3d(3 



m cosh (3 



e (mcosh/3-/x+/i if fjjO/T -r I 

1 Jo^mm^^D-^d) 



+ three terms 



4T 



T ' 



6T. 



(5.9) 



This value is twice as much as Pauli's classical value 3T obtained in the ultrarelativis- 
tic limit within the usual on-shell theory P0|, which is quite natural, since we take 
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account of particle-antiparticle pairs. The value of the specific heat is now obtained 
from (5.6), (5.9): 

(E-m)=(6-^T = ^T, (5.10) 

/ 157T 

7 =4 « 1.055. 5.11 

16 

We see that at low temperatures, 7' provides a correction of the order of 5.5% to the 
classical value 3/2. 



6 Concluding remarks 

We have discussed possible consequences of a manifestly covariant relativistic sta- 
tistical mechanics for hadronic physics. We have considered the relativistic event- 
anti-event system, for which we have obtained the relativistic mass distribution and 
calculated all characteristic thermodynamic variables. We have shown that in the 
case of small hk and zero net particle charge, the equation of state for such a system 
corresponds to the "realistic" equation of state proposed by Shuryak [4] to describe 
the behavior of hot hadronic matter. We have seen that the mass distribution ob- 
tained in this paper can serve as a framework for fitting the hadronic resonances, in 
good agreement with the phenomenological models of Hagedorn et al [8]- [10]. 

In ref. [41] we have argued that the mass potential hk is related to the degeneracy 
parameter 5 : hk ~ 8~ a , a > 0. Therefore, small fix indicates large S, i.e., high 
degeneracy. In this case the system should exhibit the collective behavior determined 
by a relatively strong interaction, mainly of attractive character. This attraction 
considerably reduces the pressure, resulting in a value for the sound velocity c 2 ~ 0.20. 

Further physical consequences of the theory of the off-shell events are discussed in 
refs. |$7|,[41]. In []37[ the statistical mechanics of the bosonic event-anti-event system 
is considered. For such a system, at some critical temperature a special type of Bose- 
Einstein condensation sets in, which provides the events making up the ensemble a 
definite mass and represents, in this way, a phase transition to an equilibrium on-shell 
sector. In [12] an adiabatic equation of state, p oc N^ 5 , is obtained for the system of 
degenerate off-shell fermions and possible implications in astrophysics are discussed. 

The various aspects of the theory of the off-shell events discussed in this paper, 
as well as related questions, are now under further investigation. 
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Appendix A 

The polylogarithm is defined by the relation 

oo 

Liy(z) = ^—, \z\<l or |z| = 1, Re v > 1. (A.l) 
k=i kv 

For \z\ > 1, the function Li v (z) is defined as the analytic continuation of this series. 
For integer n and positive z, the following relation holds (ref. [33], p. 763): 

«„ w+( _ irii „ ( i)._(M: B „gi), (A . 2) 

where are the Bernoulli polynomials (ref. [33], p. 765), 

B (z) = l, B 1 (z)=z-^ B 2 {z)=z 2 -z+ 1 -, 

B 3 (z) = z 3 -^z 2 + ^z, B A {z) = z i -2z 3 + z 2 -^, etc. 
possessing the property 



B' n (z)=nB n _ 1 (z). 



In view of (A. 2) we have 



Lii(z) — Li 1 ( - ) = — In z + ni, (A. 3) 



Li 2 (z) + Li 2 (^j = -i^ + Trilnz + y, (A.4) 
r ■ t \ t ■ (^\ ln 3 z fi, „ 7r 2 

L, s w - l.3 y = — 6~ + y y ' ^ 

/ v T / 1 \ lll 4 Z JTJ, , 7T 2 , 9 7T 4 . , , 

L. 4 («) + L. 4 ( k -J = -— + -ln 3 * + T ln 2 z + -, (A.6) 

r • / \ t ■ ( 1 \ m5 z ™ , 4 71-2 , 3 71-4 , / . j_\ 

j- . / \ r . /1\ In 6 z ixi , 5 tt 2 4 tt 4 , 2 2tt 6 /4 oX 

W,)+L. 6 y = -_ + _lnS 2 + -l^+-h^ + — .etc (A.8) 
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For integer n and negative z, we use another relation (ref. [33], p. 763): 



, i \ [«/2] 

Li n (z) + (-l) n Li n (-)=-£ c M ln B - 2 *( 
V2/ fc=i 



(A.9) 



1 7T 

n! o(n — 2)\ 



Cfe, n = (-1) 



k+1 2{2^-l)^ k B 2k 
(2k)\(n-2k)\ 



, fc = 0,l,...,[n/2], 



where B 2 k are the Bernoulli numbers, 

Bo = 1, -B 2 = -, B 4 
6 

In this way we obtain 



- — , Br = — , etc. 
30 42' 



IAi(z) 


— Lii 


Li 2 (z) 


+ Li 2 


Li 3 (z) 


- Li 3 


Li±{z) 


+ Li 4 


Li 5 (z) 


- Li 5 


Li 6 (z) 


+ Li 6 



- = -H-z), 



ln 2 (- 


z) 


2 




ln 3 (- 


z) 


6 




ln 4 (- 


z) 


24 




ln 5 (- 


z) 



120 
ln 6 (-z 



+ 



7T 



7T 



+ -ln(-z) 



H n 2 (-z) H 

12 v ; 360 

71-2 i 3 / \ 77T 4 , 

+ 36 ln( -^ ) + 360 ln( - 2) 



/T ~ , 4 7~ 2 , , 31"7T 



720 1 144 ln( ^ ) + 720 ln( ^ )4 15120 



(A.10) 
(A.ll) 

(A.12) 

(A.13) 

(A.14) 

(A.15) 



etc. 



The formulas in the main text are obtained from (A.3)-(A.8) for z = er, and from 
(A.10)-(A.15) for z = — e^, for the bosonic and the fermionic cases, respectively. 
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